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Abstract 

Even within a single genus, such as Drosophila, cases of lineage-specific adaptive evolution have been found. Therefore, the 
molecular basis of phenotypic variation must be addressed in more than one species group, in order to infer general 
patterns. In this work, we used D. americana, a species distantly-related to D. melanogaster, to perform an F2 association 
study for developmental time (DT), chill-coma recovery time (CRT), abdominal size (AS) and lifespan (LS) involving the two 
strains (H5 and W1 1) whose genomes have been previously sequenced. Significant associations were found between the 43 
large indel markers developed here and DT, AS and LS but not with CRT. Significant correlations are also found between DT 
and LS, and between AS and LS, that might be explained by variation at genes belonging to the insulin and ecdysone 
signaling pathways. Since, in this F2 association study a single marker, located close to the Ecdysone receptor {EcR) gene, 
explained as much as 32.6% of the total variation in DT, we performed a second F2 association study, to determine whether 
large differences in DT are always due to variation in this genome region. No overlapping signal was observed between the 
two F2 association studies. Overall, these results illustrate that, in D. americana, pleiotropic genes involved in the highly- 
conserved insulin and ecdysone signaling pathways are likely responsible for variation observed in ecologically relevant 
phenotypic traits, although other genes are also involved. 
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Introduction 

The complete description of the molecular basis of variation for 
quantitative traits in natural populations of different species is 
essential to understand how genetic variation for adaptive traits is 
maintained, in general, and whether the adaptive phenotypic 
variation observed within and between species is caused by the 
same variable loci [1]. In Drosophila, most of the work regarding 
this issue has been done using species from the melanogaster group 
(see for instance, [2-8]). However, other groups of species must be 
studied as well since, even within a single genus such as Drosophila, 
cases of lineage-specific adaptive evolution have been found (see 
for instance, [9-13]). Moreover, genes that have been reported as 
harboring variability that explains within species phenotypic 
variation in D. melanogaster have been found to be missing in 
distandy related Drosophila species (see for instance, [14]). 

Drosophila americana (virilis group) is becoming an important 
model for comparative studies. This species is distantly related to 
D. melanogaster since the two lineages have been diverging for about 
40 million years [15]. The large effective population size and the 
geographical distribution along the USA imply a large amount of 
nucleotide and phenotypic variation. Since D. americana is easy to 



collect across the entire range of the distribution and it is easy to 
maintain in the laboratory, a large number of strains are already 
available. Little population structure was observed for this species 
besides that caused by chromosomal inversions and a chromo- 
somal fusion, for which molecular markers are available. The 
genome of strains H5 and W 1 1 has been sequenced and assembled 
(for a detailed discussion see Fonseca et al. [16]). 

Here, we present the results of one F2 association study 
involving the two D. americana strains whose genome has been 
sequenced (H5 and Wl 1; [16]), and four phenotypic traits that are 
likely ecologically important (developmental time (DT), chill-coma 
recovery time (CRT), abdominal size (AS) and lifespan (LS)). 
Significant correlations have been found between DT and LS and 
AS and LS that can be explained by variation at genes belonging 
to the insulin and the ecdysone signaling pathways. Since, in this 
F2 association study a single marker, located close to the Ecdysone 
receptor (EcR) gene, explained as much as 32.6% of the total 
variation in DT, we performed a second F2 association study 
where the slow developing strain Wl 1 was replaced by another 
slow developing strain (W37), to determine whether large 
differences in DT are always due to variation in this genome 
region. No overlapping signal was observed between the two F2 
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Figure 1. Schematic representation of the statistics computed for the F2 association studies involving the H5<7xW11q cross for 
developmental time (DT) and life span (LS). Black bars show the difference (in percentage) between genotypic classes that show the highest 
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and lowest values for these two phenotypic traits ((highest value class - lowest value class)/lowest value class). The values next to the bars are the 
percentage of variation explained by each marker (the square of the Pearson correlation coefficient). For this calculation, we use the mean of each 
genotypic class as the expected value for that class. Significant associations are represented by *(0.05>P>0.01), **(0.01>P>0.001) and ***(P<0.001). 
Muller's elements A - E are shown in A) to E), respectively. 
doi:1 0.1 371 /journal.pone.0086690.g001 



association studies. Overall, these results illustrate the likely role of 
pleiotropic genes involved in the insulin and ecdysone pathways, 
but also of other genes, in the setting of likely ecologically relevant 
phenotypic traits. 

Materials and Methods 

Marker Development 

A set of 43 large indel markers distributed along the five major 
chromosomes of D. americana was developed based on the available 
genome sequence for strains H5 and Wll [16]. The markers 
locations, as well as the variants that are segregating in the strains 
(H5, Wll and W37) used in the F2 association studies 
(H50-xW110 and H50*xW370) are shown in Fig. SI. Primers 
and PCR amplification conditions used in this work are shown in 
Table SI. All PGR products were visualized on a UV transillu- 
minator after electrophoresis using SGTB buffer (Grisp, Portugal) 
in 2% agarose gels stained with ethidium bromide. 

F2 association Studies 

The three isofemale strains used in the F2 association 
experiments were established with flies collected at the end of 
July and beginning of August 2004 at Lake Wappapelo, Missouri 
(Wl 1 and W37) and at Lake Hurricane, Mississippi (H5) and have 
been kept in the lab at room temperature (20-22°C) in large vials 
with about 20-50 individuals. Then, these strains were transferred 
to a controlled temperature chamber at 25°C, with 12/12 hours of 
light/ dark cycles, one generation before starting the experiment. 
The F2 association experiments have been performed under the 
latter conditions. The D. americana F2 association study involving 
strains H5 and Wll is described in detail by Reis et al. [14]. 
Briefly, 131 F2 D. americana males showing extreme phenotypes, 
after excluding all individuals that show phenotypic values for any 
of the four traits (DT, CRT, AS and LS) in the second third of the 
distribution (intermediate values), were selected out of 975 
individuals. The first trait to be measured was DT. For this 
purpose, each of the 83 second generation crosses (Fl) were 
transferred to new flasks every day in order to obtain the precise 
period of time, in days, between oviposition and adult eclosion. 
The resulting F2 males were then individually collected. When F2 
males were 10 days old (young adult flies), individual CRT was 
measured, in seconds, at 25°C after four hours of cold exposure at 
0°C. Flies must be able to stand up on their legs in order to be 
considered completely recovered. Individual photographs were 
taken when individuals were 20 days old, using a stereomicroscope 
Nikon ZMS 1500 H. The resulting JPG files were saved with a 
resolution of 1600x1200 pixels. Relative AS was estimated by 
counting the number of pixels in the picture that correspond to 
this structure, using Adobe Photoshop H (Adobe, USA). The flies 
were then transferred to new vials and kept until they died, in 
order to measure LS, in days. 

For the F2 association study involving the H5 and W37 strains, 
1 30 individuals showing 14 or less days and 18 or more days of DT 
(time between oviposition and adult eclosion) were selected. These 
three strains (H5, Wl 1 and W37) were selected, since they showed 
remarkable differences regarding the studied phenotypic traits and 
they show the same chromosomal rearrangements. 



For both association studies only males were used, in order to 
avoid a possible sex effect that might exist. Genomic DNA was 
then extracted using the OJAamp DNA Mini Kit from QIAGEN 
(Izasa Portugal, Lda.) according to the manufacturer's instructions, 
and the individuals were genotyped using the 43 indel markers 
described above. 

Genotype - phenotype associations were tested using non- 
parametric tests as implemented in SPSS Statistics 19.0 (SPSS 
Inc., Chicago, Illinois) and significant values were corrected using 
the sequential Bonferroni correction for multiple testing. Using the 
same software, non-linear as well as linear regression analyses 
(including a constant) were performed between the phenotype 
observed for each genotyped individual and the expected 
phenotype (average of the phenotype obtained for each genotypic 
class), in order to estimate the amount of phenotypic variation 
explained by variation in the indel markers. 

Candidate Genes 

Gene Ontology was used as a proxy to identify D. melanogaster 
putative candidate genes for the four phenotypic traits measured in 
this study. Then, using FlyBase (http://flybase.org), we identified 
the D. virilis orthologous genes, since this species is the most 
closely-related to D. americana (approximately 4.1 My of indepen- 
dent divergence [17]) with an annotated genome, and for which 
information on orthologous genes is available. In the few cases 
where orthology information is unavailable, putative orthologous 
genes were identified using BLASTP search and the D. melanogaster 
protein sequences as query. 

The general GO terms selected were GO: 0007476 (imaginal 
disc-derived wing morphogenesis) for DT, GO: 0008340 (deter- 
mination of adult lifespan) for LS, GO: 0040014 (regulation of 
multicellular organism growth) for AS and GO: 0009409 (response 
to cold) for CRT. 

There are 318 genes included in the GO term: imaginal disc- 
derived wing morphogenesis, along the major D. melanogaster 
chromosomes, but eight of these genes [amn; CGI 67 8; CGI 27 17; 
salm; CGI 1226; SF1; smp-30; tal-lA) have no gene correspondence 
in D. virilis. Moreover, three orthologous genes (GJ1S461/ Stim; 
GJ22504/Pka-C3; GJ22515/crc) in D. virilis are present in small 
scaffolds not assigned to any of the Muller's elements making it 
hard to determine their cytological location. After excluding these 
genes there are 307 putative candidate genes for DT differences in 
D. americana. 

A total of 147 genes are associated with the GO term: 
determination of adult life span, but eight of these genes 
(CG11700; mthl3; Acp62F; mthl2; mthW; mthU; Indy-2; mthll2) have 
no gene correspondence in D. virilis. bmm and jwd genes show two 
homolog sequences each, in D. virilis and were included in the list 
of putative candidate genes. Moreover, two orthologous genes 
[GJ15405/Cdk5alpha; GJ14992/ 'chico) in D. virilis are present in 
small scaffolds not assigned to any of the Muller's elements. After 
excluding/ introducing these genes there are 139 putative candi- 
date genes for LS differences in D. americana. 

The GO term: regulation of multicellular organism growth 
comprises 41 genes but one ortholog (GJ 14992/ chico) in D. virilis is 
present in a small scaffold not assigned to any of the Muller's 
elements. After excluding this gene there are 40 putative candidate 
genes for AS differences in D. americana). 
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Figure 2. Schematic representation of the statistics computed for the F2 association studies involving the H5o>xW119 cross for 
abdominal size (AS) and chill-coma recovery time (CRT). Black bars show the difference (in percentage) between genotypic classes that show 
the highest and lowest values for these two phenotypic traits ((highest value class - lowest value class)/lowest value class). The values next to the 
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bars are the percentage of variation explained by each marker (the square of the Pearson correlation coefficient). For this calculation, we use the 
mean of each genotypic class as the expected value for that class. Significant associations are represented by *(0.05>P>0.01), **(0.01>P>0.001) and 
***(P<0.001). Muller's elements A - E are shown in A) to E), respectively. 
doi:1 0.1 371 /journal.pone.0086690.g002 



Finally, for CRT there are 1 1 genes associated with the GO 
term: response to cold but three of these genes (brv-3; brv-2; smp-30] 
have no correspondence with any D. virilis gene and therefore, only 
eight putative candidate genes could explain CRT differences in 
D. americana. 

Results 

HScyxWIlo F2 Association Study 

In the F2 association study involving strains H5 and Wll, 
significant associations were found between the large indel markers 
here reported and DT, AS and LS (Figs. 1 and 2; Tables S2-S4). 
All associations between indel markers and CRT were non- 
significant after Bonferroni correction for multiple testing (Fig. 2; 
Table S5). 

For DT, there are nine molecular markers showing significant 
associations after applying the sequential Bonferroni correction for 
multiple comparisons, one on Muller's element D and eight on 
Muller's element C (Fig. 1). The strongest signal comes from 
Muller's element C where 64 out of 307 (20.8%) putative 
candidate genes fall into this region (Figs. 1 and 2; Table S6). 
The molecular markers showing significant associations explain as 
much as 36.97% of the total variation in DT found in this cross 
(see the C5 marker; Fig. 1 and Table S2; for this marker the 
difference between the genotypic classes showing extreme values 
for DT is 4.01 days (14.83 vs 18.84 days)). When using all markers 
showing a significant association with DT after Bonferroni 
correction, and stepwise regression, as much as 52.4% (molecular 
markers C4+C6+D1) of the total phenotypic variation in DT is 
explained. 

In the case of LS, 18 associations remained significant after 
sequential Bonferroni correction, four on Muller's element D, five 
on Muller's element E and nine on Muller's element C (Fig. 1). 
Twenty six, 18 and 26 putative candidate genes out of 139 (50.4%) 
are located in these chromosomal regions (Figs. 1 and 2; Table S6). 
The molecular markers showing significant associations explain as 
much as 22. 18% of the total variation in LS found in this cross (see 
the E6 marker; Fig. 1 and Table S3; for this marker the difference 
between the genotypic classes showing extreme values for this trait 
is 30.18 days (41.89 vs 72.07 days)). When using all markers 
showing a significant association with LS after Bonferroni 
correction, and stepwise regression, as much as 50.3% (molecular 
markers D8+C7+E3) of the total phenotypic variation in LS is 
explained. 

Twenty six indel markers show significant associations with AS, 
namely six on each of the Muller's element D and E and seven on 
each of the Muller's elements A and B (Fig. 2). Thirty five putative 
candidate genes out of 40 (87.5%) are localized at the chromo- 
somal regions showing significant associations (Figs. 1 and 2; Table 
S6). As much as 18.4% of the total variation in AS present in this 
cross is explained by a single molecular marker (see the E7 marker; 
Fig. 3 and Table S4; for this marker the difference between 
genotypic classes showing extreme values for AS is 0.20 relative 
units (0.85 vs 1 .05 relative units where 1 .00 represents the mean of 
the AS measured for all the phenotyped individuals). When all the 
markers showing a significant association with AS after Bonferroni 
correction are included in the stepwise regression analysis, as 
much as 33.4% (molecular markers E5+A4) of the total phenotypic 
variation in AS is explained. 



It should be noted that some of the phenotypic traits analyzed 
here are significantly correlated (DT vs LS, Non-parametric 
Spearman's correlation = -0.447; P<0.001 and AS vs CRT; 
Non-parametric Spearman's correlation = —0.302; P<0.001). 
Correlations are also observed when the complete set of 975 
phenotyped individuals is used rather than the 131 genotyped 
individuals (DT vs LS, Non-parametric Spearman's correla- 
tion =-0.237; P<0.001; AS vs CRT; Non-parametric Spear- 
man's correlation = —0.160; P<0.001 and LS vs AS Non- 
parametric Spearman's correlation = —0.176; P<0.001). This 
could be an indication that some genes account for variation in 
more than one of the traits analyzed here. Indeed, there are five 
genes that appear as putative candidates for more than one trait 
(Table 1). One of them (EcR located on Muller's element C), is a 
putative candidate for variation in both DT and LS. The other 
four (ilp2, Sk6, InR and Nfl; the first two genes are located on 
Muller's elements D and the other two on Muller's element E) are 
putative candidates for both LS and AS. EcR (the only pleiotropic 
gene that according to our F2 association study could explain DT 
variation) is located at ~0.8 Mbp from the indel marker C8 which 
alone explains as much as 32.6% of the total phenotypic variation 
for DT. The other four genes are putative pleiotropic candidates 
genes for AS, and the indel markers closer to these genes are D4 
(—0.4 and ~ 1 .4 Mbp from ilp2 and S6k respectively) and E6 (~ 1 .5 
and ~2.4 Mbp from Nfl and InR, respectively). Using stepwise 
regression analysis, the marker on Muller's element D is excluded 
from the model and the marker E6 explains, as much as 18.7% of 
the total variation observed for AS in the F2 association study. 
Finally, all five pleiotropic genes [EcR, ilp2, Sk6, InR and Nfl) were 
implicated in LS. Using stepwise regression analysis, the molecular 
marker D4 is excluded from the model, and the indel marker C8 
along with E6 explain as much as 37.7% of the variation observed 
for LS. 

H50-XW379 F2 Association Study 

In the HSCxWl 1Q F2 association study, a single marker, 
located close to the Ecdjsone receptor (EcR) gene, explained as much 
as 32.6% of the total variation in DT. Therefore, in order to 
address whether large differences in DT are always due to 
variation in this genome region we performed a second F2 
association study where the slow developing strain Wl 1 was 
replaced by another slow developing strain (W37). This was 
feasible since 25 out of the 43 indel markers (58. 1 %) developed 
based on the H5 and Wll genomes were also segregating in the 
H50"xW37Q cross. Eight significant associations were found 
between the molecular markers and DT, namely two on each of 
the Muller's elements A, B, D and E (Fig. 3). Markers on Muller's 
element E explain as much as 16.3% of the total variation in DT 
found in this cross (see the marker E7; Fig. 3 and Table S7; for this 
marker the difference between the genotypic classes showing 
extreme values for DT is 2.79 days (14.46 vs 17.25 days)). When 
using all markers showing a significant association with DT after 
Bonferroni correction, and stepwise regression, as much as 40.5% 
(molecular markers E7+A7+E2) of the total phenotypic variation 
in DT is explained. Interestingly, none of these eight significant 
associations co-localized with those previously found in the other 
F2 association study. 
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Figure 3. Schematic representation of the statistics computed for DT on the two F2 association studies (H5ctxW119 and 

H5o>xW379). Black bars show the difference (in percentage) between genotypic classes that developed faster and slower ((fast DT class - slow DT 
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class)/slow DT class). The values next to the bars are the percentage of variation explained by each marker (the square of the Pearson correlation 
coefficient). For this calculation, we use the mean of each genotypic class as the expected value for that class. Significant associations are represented 
by *(0.05>P>0.01), **(0.01>P>0.001) and ***(p<0.001). n.a. - not applicable, the markers were not segregating in this F2 association study. Muller's 
elements A - E are shown in A) to E), respectively. 
doi:1 0.1 371 /journal.pone.0086690.g003 



Discussion 

In the H50"xW119 F2 association study presented here, flies 
were phenotyped for four traits (DT, CRT, AS and LS) and 
associations were found between markers and all traits with the 
exception of CRT. It should be noted that significant correla- 
tions were found between DT and LS and AS and LS, which 
suggests the involvement of pleiotropic genes in the setting of 
these phenotypic traits. This possibility is further supported by 
the observation of significant associations between the same 
genome region and different traits. In D. melanogaster, DT has 
been found to be correlated with other life-history traits, such as 
adult weight at eclosion [18], adult size [18-21], pre-adult 
survival [19,20] and longevity [22]. Moreover, a negative 
correlation between body size and longevity has been described 
for some mammals [23,24], although for Drosophila the relation- 
ship between these two phenotypes is less obvious (see for 
instance [25]). 

There are plenty of candidate genes reported to be associated 
with variability in DT, AS and LS in D. melanogaster (see for 
instance [26-28]). When using general GO terms, 169 D. 
americana genes could in principle explain the observed geno- 
type-phenotype associations. Although these represent 34.8% of 
the total number of putative candidate genes, they are still too 
many to study in detail. Five genes (EcR, ilp2, S6k, InR and JVfl) 
can in principle account for variation in more than one trait. InR, 
ilp2 and S6k are members of the insulin-signaling pathway that 
along with Mfl have been implicated in the determination of body 
size [28-30] and adult life span [31-34]. EcR has been shown to 
have a role in development during metamorphosis [35] and it is 
involved in the regulation of longevity [36]. These five genes 
might explain as much as 18.7, 32.6 and 37.7% of variation in 
AS, DT and LS, respectively in this particular cross. This is 
respectively 56.0%, 62.2% and 75.0% of the variation that is 
explained when using all markers that show a significant 
association with these traits after Bonferroni correction. There- 
fore, in Drosophila, genus-wide, a fraction of the variation in likely 
ecologically relevant traits (such as DT, AS and LS) might be due 
to pleiotropic genes. 

In the H50"xWl 19 F2 association study a significant association 
was obtained for almost the entire Muller's element C. In D. 
melanogaster, Mensch et al. [37] found natural allelic variants that 
contribute to variation in DT for the genes mastermind; invected, 



cricklet and CG14591, that are all located on Muller's element C. 
Therefore the result we get for D. americana could be taken as 
evidence for the presence of several genes influencing DT in this 
Muller's element. Nevertheless, we cannot rule out the possibility 
of the involvement of a single gene of major effect associated with 
a low recombination rate along Muller's element C as the 
explanation for our result. One such gene could be EcR that is 
located close to a marker that explains as much as 32.6% of the 
total variation in DT. In order to determine whether large 
differences in DT are always due to variation in this genome 
region we performed a second F2 association study where the slow 
developing strain Wl 1 was replaced by another slow developing 
strain (W37). None of the statistically significant genotype- 
phenotype associations found in this second F2 association study 
co-localized with those obtained in the first one. Therefore, EcR 
cannot explain DT differences in the second study. One possibility 
is that the experimental design of the first study (selecting only 
individuals that show extreme phenotypes simultaneously for the 
four traits surveyed) favors the identification of this pleiotropic 
gene. Indeed, the correlation between DT and LS is higher when 
the set of 131 individuals showing extreme phenotypes (Non- 
parametric Spearman's correlation = —0.447; P<0.001) is com- 
pared with the entire dataset (975 individuals; Non-parametric 
Spearman's correlation = -0.237; P<0.001). Huang et al. [38] 
have reported that epistasis is an extremely important factor at 
defining variation for quantitative traits. Thus, another possible 
explanation is that the use of a different strain showing different 
alleles may have disrupted the interactions between EcR and other 
genes associated with developmental time differences. Other 
epistatic interactions could have been created in the second study 
resulting in different regions showing statistically significant 
associations. Therefore, at present, it seems likely that, in D. 
americana, variation in DT is due to the involvement of the 
pleiotropic gene EcR, but also of multiple other genes located 
across the entire genome. Compatible with this view, in D. 
melanogaster, about 65% of the 178 mutant lines analyzed by 
Mensch et al. [26] showed altered DT. 

Despite the large number of candidate genes for DT, there is 
one region in the H50"xWl 19 F2 association study that shows a 
significant association with this trait and for which there is no 
candidate gene (see marker Dl; Fig. 3 and Table S6). This may 
be due to incomplete knowledge on genes that affect DT 



Table 1. List of putative candidate genes that could account for variation in more than one phenotypic trait (pleiotropic genes). 





Gene 


Location (£>. virilis) 


In between 
molecular markers 


DT 


LS 


AS 


EcR (C) 


scaffoldj 2875:4,666,484.4,709,770 


C7-C8 


0.000*** 


0.000*** 




Hp2 (C) 


scaffoldj 3049:8,664,534.8,665,005 


D4-D5 




0.007** 


0.000*** 


S6k (D) 


scaffoldj 3049:9,771 ,91 1 .9,782,583 


D4-D5 




0.007** 


0.000*** 


InR (D) 


scaffoldj 3047:6,346,220.6,353,846 


E5-E6 




0.000*** 


0.000*** 


Nfl (E) 


scaffold J 3047:5,385,397.5,395,947 


E5-E6 




0.000*** 


0.000*** 



The Muller's element is shown within brackets after the gene name. Significant associations are denoted by **(0.001 <P<0.01 ) and ***(P<0.001 ). It should be noted that 
the reported associations are with the markers (see Tables S2-S4) and not with the genes. 
doi:1 0.1 371 /joumal.pone.0086690.t001 
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variation in Drosophila, but we cannot rule out the possibility that 

D. americana lineage specific genes are responsible for such an 
association. Such lineage-specific genes can only be identified 
using an approach that does not rely on candidate genes 
identified in other species. 

It could be argued that the latitudinal and altitudinal gradients 
for DT that are found in D. melanogaster [39-41] and in D. buzzatii 
[42], and that are likely also present in D. americana (see [16]) could 
imply the presence of at least some common variants influencing 
DT in natural populations. Therefore, it is still conceivable that 
common variants will be found in D. americana when performing 
F2 association studies involving individuals from the extremes of 
the distribution. Nevertheless, performing such association studies 
may be difficult due to the suppression of recombination effect 
caused by the presence of the XI 4 fusion (Muller's elements A/B) 
and six inversions (Xc, 2b, 4ab, 5a and 5b; on Muller's elements A, 

E, B, C and C, respectively [43]), that are found at very different 
frequencies in populations from the extremes of the distribution. 

Overall, our results illustrate not only the role of pleiotropic 
genes belonging to the insulin and ecdysone pathways, but also of 
other genes in the setting of likely ecologically relevant phenotypic 
traits. The set of 43 cheap and easy to type large indel markers 
distributed along the five major chromosomes of D. americana that 
was developed here is a tool that will ease the identification of the 
causative genes. 

Supporting Information 

Figure Si Schematic representation of the indel mark- 
ers developed for all five large D. americana chromo- 
somal arms (Muller's elements A - E are shown in A) to 
E), respectively). The order is the same in every panel: from left 
to the right is shown the 1 Kb DNA ladder (GeneRuler™ 1 Kb 
DNA ladder from Thermo Scientific, USA), the H5, Wl 1 and 
W37 strains. Some of the variants found in the H5 individuals 
involved in the H50" x Wl 1Q cross were not segregating in the 

References 

1. Mackay TFC (2001) Quantitative trait loci in Drosophila. Nat Rev Genet 2: 1 1- 
20. 

2. Catania F, Kaucr MO, Dabom PJ, Yen JL, Ffrench-Constant RH, et al. (2004) 
World-wide survey of an Accord insertion and its association with DDT resistance 
in Drosophila melanogaster. Mol Ecol 13: 2491-2504. 

3. Schlenke TA, Begun DJ (2004) Strong selective sweep associated with a 
transposon insertion in Drosophila simulans. Proc Natl Acad Sci U S A 101: 1626- 
1631. 

4. Aminetzach YT (2005) Pesticide resistance via transposition-mediated adaptive 
gene truncation in Drosophila. Science 309: 764—767. 

5. Pool JE, Aquadro CF (2007) The genetic basis of adaptive pigmentation 
variation in Drosophila rnclcuwgaslcr. Mol Ecol 16: 2844—2851. 

6. Bono JM, Matzkin LM, Castrezana S, Markow TA (2008) Molecular evolution 
and population genetics of two Drosophila mcllkri cytochrome P450 genes involved 
in host plant utilization. Mol Ecol 17: 321 1-3221. 

7. Matzkin LM (2008) The molecular basis of host adaptation in cactophilic 
Drosophila: molecular evolution of a Glutathione 5-transferase gene (GstDl) in 
Drosophila mojavensis. Genetics 178: 1073-1083. 

8. Fry SN, Rohrseitz N, Straw AD, Dickinson MH (2008) TrackFly: Virtual reality 
for a behavioral system analysis in free-flying fruit flies. J Ncurosci Methods 171: 
110-117. 

9. Barbash DA, Awadalla P, Tarone AM (2004) Functional divergence caused by 
ancient positive selection of a Drosophila hybrid incompatibility locus. PLoS Biol 
2: cl42. 

10. Presgraves DC, Stcphan W (2007) Pervasive adaptive evolution among 
interaetors of the Drosophila hybrid inviability gene, Nup96. Mol Biol Evol 24: 
306-314. 

1 1 . Bachtrog D (2008) Similar rates of protein adaptation in Drosophila rniranda and 
D. melanogaster, two species with different current effective population sizes. BMC 
Evol Biol 8: 334. 

12. Llopart A, Comcron JM (2008) Recurrent events of positive selection in 
independent Drosophila lineages at the spermatogenesis gene roughex. Genetics 
179: 1009-1020. 



H5C x W37Q cross (markers B2, B5, C6, D8, E3) and therefore 
these molecular markers could not be used in the second study. 
(TIF) 

Table SI Primers used for the 43 equidistant molecular 
markers distributed along the D. americana genome 
based on the two completely sequenced D. americana 
strains (H5 and Wll; Fig. SI). 

(PDF) 

Table S2 Results of the F2 H50"xWll9 association study 
for developmental time. 

(XLS) 

Table S3 Results of the F2 H50-XW119 association study 
for life span. 

(XLS) 

Table S4 Results of the F2 H50"xWl 19 association study 
for abdominal size. 

(XLS) 

Table S5 Results of the F2 H50-XW119 association study 
for chill coma recovery time. 

(XLS) 

Table S6 List of candidate genes located near the 
chromosomal regions showing statistically significant 
associations after Bonferroni correction. 

(XLSX) 

Table S7 Results of the F2 H5o-xW379 association study 
for developmental time. 

(XLSX) 

Author Contributions 

Conceived and designed the experiments: MR CPV JV. Performed the 
experiments: MR IP HR BA CPV JV. Analyzed the data: MR CPV JV. 
Contributed reagents/matcrials/analysis tools: MR CPV JV. Wrote the 
paper: MR IP HR BA CPV JV. 



13. Morales-Hojas R, Vicira CP, Reis M, ViciraJ (2009) Comparative analysis of 
five immunity-related genes reveals different levels of adaptive evolution in the 
virilis and melanogaster groups of Drosophila. Heredity 102: 573—578. 

14. Reis M, Vieira CP, Morales-Hojas R, Aguiar B, Roeha H, et al. (2011) A 
comparative study of the short term cold resistance response in distantly related 
Drosophila species: The role of regucalcin and Frost. PLoS ONE 6: c25520. 

15. Morales-Hojas R, Vieira J (2012) Phylogenetic patterns of geographical and 
ecological diversification in the subgenus Drosophila. PLoS ONE 7: c49552. 

16. Fonseca NA, Morales-Hojas R, Reis M, Rocha H, Vieira CP, et al. (2013) 
Drosophila americana as a model species for comparative studies on the molecular 
basis of phenotypic variation. Cenome biology and evolution 5: 661-679. 

17. Morales-Hojas R, Reis M, Vieira CP, ViciraJ (201 1) Resolving the phylogenetic 
relationships and evolutionary history of the Drosophila virilis group using 
multilocus data. Mol Phylogen Evol 60: 249-258. 

18. Nunney L (1996) The response to selection for fast larval development in 
Drosophila melanogaster and its effect on adult weight: an example of a fitness trade- 
off Evolution 50: 1193. 

19. Chippindale AK, Alipaz JA, Chen H-W, Rose MR (1997) Experimental 
evolution of accelerated development in Drosophila. 1. Developmental speed and 
larval survival. Evolution 51: 1536. 

20. Prasad NO, Shakarad M, Anitha D, Rajamani M, Joshi A (2001) Correlated 
responses to selection for faster development and early reproduction in 
Drosophila: the evolution of larval traits. Evolution 55: 1363—1372. 

21. Zwaan B, Bijlsma R, Hoekstra RE (1995) Artificial selection for developmental 
time in Drosophila melanogaster in relation to the evolution of aging: direct and 
correlated responses. Evolution 49: 635. 

22. Chippindale AK, Alipaz JA, Rose MR (2004) Experimental evolution of 
accelerated development in Drosophila. 2. Adult fitness and the fast development 
syndrome. Methuselah Hies: a case study in the evolution of aging World 
Scientific Press, Singapore: 413—435. 

23. EigenmannJE, Patterson DE, Eroesch ER (1984) Body size parallels insulin-like 
growth factor I levels but not growth hormone secretory capacity. Acta 
Endocrinol (Copenh) 106: 448-453. 



PLOS ONE | www.plosone.org 



8 



January 2014 | Volume 9 | Issue 1 | e86690 



D. americono Pleiotropic Genes 



24. Miller RA, Chrisp C, Atchley W (2000) Differential longevity in mouse stocks 
selected for early life growth trajectory. J Gerontol a-Biol 55: B455— B461. 

25. Clancy DJ, Gems D, Harshman LG, Oldham S, Stacker H, et al. (2001) 
Extension of life-span by loss of CHICO, a Drosophila insulin receptor substrate 
protein. Science 292: 104-106. 

26. Mensch J, Lavagnino N, Garrcira V, Massaldi A, Hasson E, et al. (2008) 
Identifying candidate genes affecting developmental time in Drosophila melano- 
gaster. pervasive plciotropy and gene-by-environment interaction. BMC 
Developmental Biology 8: 78. 

27. Paaby AB, Schmidt PS (2009) Dissecting the genetics of longevity in Drosophila 
melanogaster. Fly 3: 29-38. 

28. Oldham S, Bohni R, Stocker H, Brogiolo W, Hafen E (2000) Genetic control of 
size in Drosophila. Philos Trans R Soc Lond B Biol Sci 355: 945-952. 

29. Zhang H, Liu J, Li CR, Momen B, Kohanski RA, et al. (2009) Deletion of 
Drosophila insulin-like peptides causes growth defects and metabolic abnormal- 
ities. Proc Natl Acad Sci U S A 106: 19617-19622. 

30. Walker JA, Tchoudakova AV, McKcnncy PT, Brill S, Wu D, et al. (2006) 
Reduced growth of Drosophila neurofibromatosis 1 mutants reflects a non-cell- 
autonomous requirement for GTPase-Activating Protein activity in larval 
neurons. Genes Dev 20: 3311-3323. 

31. Gronkc S, Clarke D-F, Broughton S, Andrews TD, Partridge L (2010) 
Molecular evolution and functional characterization of Drosophila insulin-like 
peptides. PLoS Genet 6: el000857. 

32. Ikeya T, Broughton S, Alic N, Grandison R, Partridge L (2009) The 
endosymbiont Wolbachia increases insulin /IGF- like signalling in Drosophila. 
Proc Biol Sci 276: 3799-3807. 

33. Kapahi P, Zid BM, Harper T, Koslover D, Sapin V, et al. (2004) Regulation of 
lifespan in Drosophila by modulation of genes in the TOR signaling pathway. 
Curr Biol 14: 885-890.' 



34. TongJJ, Schriner SE, McCleary D, Day BJ, Wallace DC (2007) Life extension 
through ncurofibromin mitochondrial regulation and antioxidant therapy for 
neurofibromatosis- 1 in Drosophila melanogaster. Nat Genet 39: 476-485. 

35. D'Avino PP, Thummcl CS (2000) The ecdysone regulatory pathway controls 
wing morphogenesis and integrin expression during Drosophila metamorphosis. 
Dev Biol 220: 211-224. 

36. Simon AF, Shih C, Mack A, Benzer S (2003) Steroid control of longevity in 
Drosophila melanogaster. Science 299: 1407-1410. 

37. Mensch J, Garrcira V, Lavagnino N, Gocnaga J, Folguera G, et al. (2010) Stage- 
specific effects of candidate hetcrochronic genes on variation in developmental 
time along an altitudinal clinc of Drosophila melanogaster. PLoS ONE 5: el 1229. 

38. Huang W, Richards S, Carbone MA, Zhu D, Anholt RRH, ct al. (2012) 
Epistasis dominates the genetic architecture of Drosophila quantitative traits. 
Proceedings of the National Academy of Sciences 109: 15553—15559. 

39. James AC, Partridge L (1995) Thermal evolution of rate of larval development in 
Drosophila melanogaster in laboratory and field populations. J Evol Biol 8: 315—330. 

40. Van't Land J, Van Putten P, Zwaan Kamping, Delden WV (1999) Latitudinal 
variation in wild populations of Drosophila melanogaster: heritabilities and reaction 
norms. J Evol Biol 12: 222-232. 

41. Folguera G, Ceballos S, Spezzi L, Fanara JJ, Hasson E (2008) Clinal variation in 
developmental time and viability, and the response to thermal treatments in two 
species of Drosophila. Biol J Linn Soc Lond 95: 233—245. 

42. Sambucetti P, Loeschcke V, Norry FM (2006) Developmental time and size- 
related traits in Drosophila Imzzatii along an altitudinal gradient from Argentina. 
Hereditas 143: 77-83. 

43. Hsu TC (1952) Chromosomal variation and evolution in the virilis group of 
Drosophila. University of Texas Publications 5204: 35-72. 



PLOS ONE | www.plosone.org 



9 



January 2014 | Volume 9 | Issue 1 | e86690 



